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Abstract 

We investigate cross-correlations in the tunneling currents through two parallel quantum dots 
coupled to independent electrodes and gates and interacting via an inter-dot Coulomb interaction. 
The correlations reveal additional information, beyond what can be learned from the current or 
conductance, about the dynamics of transport processes of the system. We find qualitatively 
different scenarios for the dependence of the cross-correlations on the two gate voltages. Reducing 
the temperature below the inter-dot Coulomb interaction, regions of a given sign change from 
spherical shapes to angular L-shapes or stripes. Similar shapes have been observed in recent 
experiments. 
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Introduction - The analysis of shot noise and current correlations in quantum dots 
reveals important information about the dynamics of the transportprocesses in the system, 
more than what can be obtained from the current or conductance, jl] Due to the presence of 
extrinsic noise sources (e.g. l//-noise) experiments are difficult to perform, but very recently 



two groups succeeded in measuring the so-ca 



flowing through a quantum dot system. 

Baa 

Cross-correlations are also important when 



led cross-correlations between two currents 



a quantum dot setup is used as a beam-splitter 

In this letter we study the cross-correlations for a model of the double dot system studied 
by McClure et al. in Ref. 2). Two parallel quantum dots (top and bottom), interacting via 
an inter-dot Coulomb interaction, are coupled to four independent source or drain electrodes 
and two independent gates. The cross-correlations between the currents through the two 
dots, plotted as a function of two gate voltages at fixed transport bias, form regions of 
positive or negative sign. Depending on the strength of the inter-dot Coulomb interaction, 
the temperature, and the transport bias we observe various shapes of these regions, namely 
spherical shapes, L-shapes, or stripes. In particular, L-shapes and stripes develop as the 
temperature is lowered, confirming what may have been observed in experiments. [7] We 
find that the observed structure of the cross-correlations is due to the combination of various 
transport processes between the charge states that can be reached at the given transport 
bias. 

Model - We consider two parallel single-level quantum dots (top and bottom (i = 
t,b)) with level energies e«, onsite Coulomb repulsion U and inter-dot Coulomb interaction 
U nn .The double dot Hamiltonian reads 



icr i <ij>aa' 



Since tunneling between the dots is suppressed the energies of the double dot eigenstates are 
fully determined by the dot occupation number n ia = c\ a c ia . The dots are coupled via the 
usual terms H^ r = J^ika yi a \(rr c ™ + ^- c - ) ^° ^ our independent source and drain electrodes 
(r = l..,4, see left panel of Fig. [[J. The electrodes (with operators dkar) are described by 
non- interacting electrons with density of states p r and chemical potential fi r . The tunneling 
amplitudes t\ are chosen independent of electron spin a and lead momentum k. 

The gate voltages Vo,t and Vc,b shift the energies of the double dot eigenstates by eVc,t • 
(n t + a nb) + eVa,b • ( n b + ol n t) where n- x = J2 a n ™ are the occupation numbers of the dots 
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Jl. 



FIG. 1: ie/i panel: sketch of the system. Two parallel quantum dots are coupled to four electrodes, 
interacting via the inter-dot Coulomb repulsion U nn . (Color online) Right panel: sketch of two 
competing/supporting processes which lead to a different sign of the cross-correlations. The dashed 
(solid) arrow indicates transport over the bottom (top) dot. 

i = t,b. The factor a is a measure for the unavoidable (in experiment) cross talk between 
gate voltages, which is of the order a — 0.1. [2J To compare to these experiments we choose 
the two left electrodes (and similar the two right ones) to be on the same potential /i^, 
and apply the transport bias V tran s = H>l — (J>r symmetrically. 

The current operators are defined by I r = —i(e/h) ^2 ik(7 (ti^k^Cia — h.c/j. Since there 
is no tunneling between the dots the currents over the top (bottom) channel are conserved. 
Therefore, we consider symmetric combinations over the corresponding source and drain 
electrodes, e.g. for current through the top dot I t = (Ji— 1 2 )/2 and = (J t ). For each current 
the (zero-frequency) shot noise power spectrums is given by Si = J^° oo dT(SIi(r)6Ii(0) + 
8Ii(0)8Ii(r)) with 5Ii(r) = Ii(r) — (ij). The cross-correlations S t b between top and bottom 
currents are defined as 

/oo 
dr(5i t (r)5I h (0) + 5i t (0)5i h (r)). (1) 
-oo 

Expressions for the shot noise and and cross correlations in the sequential tunneling limit 
were first derived by Korotkov (Ref. ^ and references therein) within a master rate equation 

n n 

approach for single electron transistors, and by Cottet et al. [5|, [6( for several quantum 
dot models. We can analyze the noise and cross-correlations in the frame of the real-time 



diagrammatic expansion lfj] in the coupling constants T r = 2ii\t r \ 2 p r to the leads. [8|, 111] 
This is useful if higher order terms such as co-tunneling are investigated. In the present 
problem the main features are found in lowest order, i.e., for sequential tunneling already. 
To describe this situation we can use Fermi's golden rule to obtain the transition rates W xa > 
between eigenstates \x) and |%') of the quantum dot system. Considering the processes 



involving the electrode r we have 



£Wl4lx> 



+ 



[1 - fr(-E x , iX )} 



(2) 



for x' 7^ X, together with W^ x = - E x ^x W x',r Here ^x',x = E x' ~ E x is tne energy 
difference between the states involved, f(x) = f/(exp (x/T) + f) is the Fermi function, and 
f r (x) = f(x — fi r ). The total rate is the sum over all electrodes, W = W r (bold face 
symbols indicate matrix/vector notation). From the transition rates we obtain the stationary 
occupation probabilities p st from the stationary master equation, Wp sl = 0. To compute the 
current we also need to account for the direction of tunneling process to/from the quantum 
dots and the considered lead. The corresponding current rate, e.g., for the top current is 
given by W^ x , = (W^, -W^ %I ){Q{N X - N x >) - Q(N X > - N x )), where N x is the total number 
of electrons in the state |x). The currents can then be expressed as = 7^e T W /t > b p st , 
where for the efficiency of the notation we introduced e — (1, . . . , 1). 8[] 

In sequential tunneling, the expression for the shot noise (auto-correlation of a given cur- 
rent) consists of two terms. One comes from the "self-correlation" of the two current opera- 
tors, the other is due to the time evolution of the quantum dot system (in the presence of the 
electrodes) in the time between the action of current operators. For the cross-correlations, 

9I Therefore, the cross-correlations are 



the first term vanishes in sequential tunneling, 
concisely written as 

p 2 

S th = -e T W 7t PW /b p st . (3) 

n 

Here, the conditional probability P accounts for the time evolution mentioned above. [9(] It 
is determined via WP = p st (g) e T — 1, which can be derived from the Dyson equation for 
the full time propagator of the quantum dot system. 8| The more complicated expressions 
for the cross-correlations in higher order perturbation theory are a straightforward general- 



ization 



12| of the shot noise expressions, though care is necessary to correctly account for 



the rates due to the various electrodes. 

Results - Transport in quantum dot systems with a large charging energy U and weak 
coupling r r is dominated by sequential processes, in which a quantum dot is charged by a 
single excess electron at a time. The transport voltage needs to overcome the sequential 
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FIG. 2: (Color online) Left panel: Stability diagram. Zero bias conductance vs. both gate voltages 
Vc,t and Vc,b- The dashed lines separate the ground states with different occupation number pairs. 
Between the four degeneracy vertices the ground state occupation is (1,1). Right panel: Cross- 
correlations vs. both gate voltages at a transport voltage Vt ra ns = O.lmV. Red (bright) regions 
indicate positive and blue (dark) regions negative cross-correlations. 

threshold for the single-electron charging, otherwise the transport is suppressed (Coulomb 
blockade). For the coupled dot system considered the sequential threshold of each dot 
depends on the gate voltages, and because of the inter-dot Coulomb repulsion U nn on the 
occupation of the other dot. The interactions between the dots may result in currents which 
mutually support or suppress each other, leading to positive or negative cross-correlations 

S'tb- 

For zero gate voltages, we consider dots with identical level energies e t = = —0.2 
and intra-dot Coulomb repulsion U = 0.2 (all energies in meV). The inter-dot repulsion 
U nn = 0.02 is on the order of the temperature k-^T = 0.024 for Fig. [2j The dot-electrode 
couplings are chosen equal, T r = 0.003 = k B T/8. This defines the line width of a dot level 
as T = 2T r = 0.006. 

The stability diagram of the system, the differential conductance dl/dV (in units of 
Te 2 /h) of the total current I = I t + lb at zero transport voltage vs. both gate voltages, is 
shown in the left panel of Fig. [2j The lines of non-zero conductance separate the regions 
with ground states corresponding to pairs of occupation numbers ((n t ), (n^)), and form a 
(slightly) tilted square lattice (honeycomb lattice, if inter-dot tunneling would be included). 
The intersection points (vertices) of the lines correspond the situation when ground states of 
four occupation number pairs are degenerate. As for our single-level model the occupation 
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is limited to at most two electrons per dot, we have nine possible occupation number pairs 
(with corresponding spin degeneracy) Between the four degeneracy vertices, the ground state 
has occupation ((n t ), (nt,)) = (l, 1), i.e. each dot is occupied by one electron at zero bias. 

We now apply a fixed transport voltage V tTans = O.lmV over each dot and plot the cross- 
correlations (also in units of Te 2 /h) of the system vs. both gate voltages, see right panel of 
Fig. [2j We find sixteen sectors with positive (red) and negative correlations (blue). The 
cross-correlations form a checkerboard pattern, where four sectors form a four-leaved clover 
shape around a degeneracy vertex. Note that the sign of the cross-correlations is not simply 
a ground state property, as there are e.g. four spherical regions of alternating sign in the area 
with ground state (1, 1). Such a clover structure around degeneracy vertices was observed 
in experiments |2| with two parallel multi-level quantum dots. The different signs were 
explained within a minimal model of two competing/ supporting processes. 

Let us consider the clover structure at the bottom left degeneracy vertex of Fig. [2] (choos- 
ing another vertex would simply need a change of occupation numbers) . In the area around 
zero gate voltages, Vo,t ~ Vc,b ~ 0, the ground state is ((n t ), (rib)) =(1, 1) and there are two 
nearly degenerate excited states (2, 1) and (1, 2) that can be reached at the given transport 
bias. The two transport processes (1,1) — > (1 ; 2) and (1,1) — > (2,1) correspond to cur- 
rents over different dots, see right panel of Fig. [TJ lower part. Both processes have the same 
probability but only one is possible at a given time, as the realization of one process preempts 
the other. Therefore, the two currents are competing and the resulting cross-correlations 
are negative. 

We now decrease the bottom gate voltage so that we reach the region of positive cross- 
correlations. For gate voltages, Vo,t ~ O.OmV, Vc,b ~ —O.lmV, the ground state is now 
(1,2), the closest excited state is (1,1). Transport via the lower dot, (1,2) — > (1,1), now 
allows transport over the upper dot, via (1, 1) — > (2, 1), which is still within reach at the 
given transport bias, see right panel of Fig. [IJ upper part. Therefore, current over the lower 
dot supports the transport over the upper dot, and the cross-correlations are positive. 

We ignored so far that because of the relatively small inter-dot interaction U nn the state 
(2, 2) can also be reached from (1, 2) (or (2, 1)) in nearly the same gate voltage region than 
the states (1,2) ((2,1)) can be reached from (1,1). Therefore, the processes involving the 
state (2, 2) are also involved. However, the above reasoning still holds if the temperature T 
is larger than U nn , which was the case so far. 
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FIG. 3: (Color online) Left panel: Sketch of the processes relevant for the lower left negative 
L — shape. Different combinations of positive and negative processes result in the L — shape of 
the cross-correlations. Right panel: Cross-correlations vs. both gate voltages Vc,t and Vc.b at a 
temperature A;bT = 0.008. The cross-correlations develop into L — shapes inside the central region 
with ground state (1, 1). The dashed lines separate the regions of different ground state occupation. 



This changes dramatically when we decrease the temperature, such that U nn = 0.02 > 
ksT = 0.008 (we also lower the lead couplings to T r = 0.001). Now we observe a very 
different structure of the cross-correlations, as seen in the right panel of Fig. [3j The 
symmetric clover shape around a degeneracy vertex changes to L-shapes in the region with 
the ground state (1, 1), and stripes/bubbles for the other regions. 

The L-shape around zero gate voltages (ground state (1, 1)) is due to a break up of the 
previously spherical region with negative cross-correlations into four regions, which can be 
distinguished by the number of ways (if any) the state (2, 2) can be reached at the given bias. 
As discussed above, around zero gate voltages the system has processes (1, 1) — > (1, 2) and 
(1, 1) — > (2, 1) available, which leads to negative cross-correlations. On the other hand, 
if the system can also reach the state (2,2) from (2,1) and (1,2), this would bring about 
positive cross-correlations. The value of the cross-correlations is a now a question of how 
much the "positive" processes can compensate the negative cross-correlations of the processes 
related to the ground state (1, 1). 

If both gate voltages are positive, the state (2, 2) can not be reached at all, and we are 
back to the earlier explanation of (maximum) negative cross-correlations. If we decrease 
one of the gate voltages (but keep the other fixed) at some point the transport voltage 



becomes sufficient so that the system can make one of the processes (1, 2) 



:2,2) or 
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FIG. 4: (Color online) Left panel: Differential conductance vs. both gate voltages Vet an d Vc,b 
at a fixed bias Vt ran s = O.lmV and k^T = 0.008. The conductance is enhanced in the areas 
with maximum positive cross-correlations. Right panel: Cross-correlations at a small bias 14 ran s = 
0.05mV vs. both gate voltages. The cross-correlations regions show clover structure and are 
non-zero only near the degeneracy vertices. 

(2, 1) — ► (2,2), depending on which gate voltage was lowered. This additional "positive" 
process leads to less negative cross-correlations in the "legs" of the L-shape. 

If we decrease both gate voltages in the same way from zero gate voltage, going along the 
diagonal in Fig. [31 we reach an area where the cross-correlations are practically vanishing. 
Here, the state (2, 2) can be reached from both states (1, 2) and (2, 1) and, in combination, 
the positive and negative processes effectively cancel for the cross-correlations. Nevertheless, 
a finite current flows over both dots. 

The L-shapes with positive cross-correlations are explained by a similar combinations of 
processes, though in reverse order, such that the initially positive cross-correlations are more 
and more compensated by negative processes as the gate voltages are changed. That we 
observe no L-shapes for the regions not associated with the (1, 1) ground state is partly due 
to the limitation of our model to single-level quantum dots. Multi-level dots with a finite 
level spacing could induce even more structure than presently observed. On the other hand, 
the present shapes are clearly a consequence of the four charge states around the degeneracy 
vertices. This structure should be fairly robust given a sufficiently large charging energy U. 

As the transport bias dictates which states are available for sequential transport its value 
will strongly influence the observed shape of the cross-correlations. In the right panel of 
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Fig. H] we display the cross-correlations vs. both gate voltages for a reduced bias voltage 
of Vtrans = 0.05mV. In this plot the L-shapes vanish and only clearly separated clover 
structures at the degeneracy points are visible, similar to the cross-correlations for higher 
temperature. At larger transport bias (not shown) the regions with finite cross-correlations 
grow until they start to "overlap" . 

The left panel of Fig. @] depicts the differential conductance dl/dV of the total current 
vs. both gate voltages for the fixed bias voltage Vt ra ns = O.lmV. The conductance has some 
similarity to the cross-correlations. Transport over each dot results in four vertical and 
horizontal lines of high conductance (orange-yellow), leading to sixteen intersection points. 
These intersection points are at same gate voltages than the points of maximum amplitude 
of the cross-correlations of Fig. [3J In addition, intersection points with maximum positive 
cross-correlations correspond to maximum conductance, whereas points with negative cross- 
correlations show lower conductance. However, there are lines with high conductance where 
the cross-correlations essentially vanish. There, the transport is carried only over one dot, 
while the other dot is in the Coulomb blockade. 

Summary. - In summary, we found that the cross-correlations of currents through two 
parallel quantum dots are strongly influenced by the inter-dot Coulomb interactions, the 
temperature and the transport bias. Processes favoring different signs of cross-correlations 
lead to various shapes, similar to experimental findings. If the inter-dot Coulomb interaction 
is much larger than the temperature, the shape changes with increasing transport voltage 
from a four leaved clover around the charge degeneracy points, to L-shapes and stripes, 
before the regions start to overlap at large bias. In contrast to the stability diagram, the 
cross-correlations are not determined by the ground state alone, but depend sensitively on 
the excited states and the sequence of processes by which these states are reached. 
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